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Abstract 



O 

. We prove that a uniformized variant of both the Rosenthal walk [10] and the Kac 

! random walk [6] on SO{n) mixes in O(n^) steps in total variation distance. The 

' proof also extends easily to Rosenthal walk with fixed angle 6 ^ tt. To the best of 

our knowledge, this is the first polynomial time bound for both walks. The techniques 
\ employed are mainly from representation theory of 50 (n). But a crucial new ingredient 

is the interpretation of the Fourier coefficients of the character ratio as counting the 
■ number of particle cascade paths arising from the classical branching rules. 



Consider the random walk on SO{N) defined by Xq = id, Xj = Xj^iSj, where Sj are 
lid random elements in SO{N) uniformly supported on the conjugacy class of -R(l, 2; 6). In 
other words, each step we take a random rotation of angle 6 along a uniformly chosen two- 
plane. The goal is to bound the total variation mixing time to its stationary distribution. 



Q\ I namely the Haar measure. 

For two probability measures on the same space fl, recall the following consequence of 
the Cauchy-Schwarz inequality: 

II ||2 _ II 1 ||2 

11/^ - '^IItv - 11^ - ^Wc^du) 

<\& 1l|2 

5^ ! Here we are assuming the density ^ exists, but even when it doesn't, the first equality still 
holds if we restrict the integral to where ^ < oo. However, the C? norm might be infinite. 

Now we specialize to f2 = SO{N) and u = U, the Haar measure. Then by Plancherel's 
theorem, the right hand side equals (see [5] Chapter 3) 

J2 Tr(/:^)(p)(/r^)(pf = 5^Tr/i(p)/I(pf. 

To see the last equality, simply observe that for p = 1, and any probability measure vr, 
7r(p) = 1, hence (p — U){p) = 0; for any other p, U{p) = since p{g)U{p) = U{p) for any 



g G SO{N), by left-invariance of U. 
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Furthermore, we have the convolution identity 'ji*l){p) = Ji{p)u{p). Now let S'/s be 
distributed according to /x. Then since /i is constant on conjugacy classes, by Schur's lemma, 
we know that is always a constant times the identity matrix. Therefore 

Tr/?*(p)/?(p)^ = dl{c,ld,f\ 

where Cp = Tr /i(p). 

Now we further restrict to the case of = 2r;, + 1. Then the irreducible representations 
of SO{N) are indexed by n-tuples (0 < ai < 02 < . . . < a„) [7]. We will denote da and 
Ca = Ca{d) the dimension and the character value evaluated at i?(l,2;^) of the irreducible 
representation labeled by the n-tuple a. Note that Rosenthal [10], [1] used the following 
slightly different labeling convention (ai + 1/2, 02 + 3/2, . . . , a„ + n — 1/2). He gave the 
following explicit formulae for da and Ca by taking suitable limits in the Weyl character 
formula: 

4 = Yi^j (2n-l)! n("'? + g- V2) n [K + r- 1/2)^ -(a, + .-1/2)2] 

^ q=l l<s<r<n 

{2n- 1)! ^ sin(aj + j - 1/2)6 

'"^''^ - (2sin(^/2))2-i 1^ (a, + J - 1/2) n.^,(«r + r - 1/2)^ - (a, + j - 1/2)^ 

By restricting to the caononical copy of SO {2) in SO{N), it decomposes into a non- 
negative linear combination of irreducible representations of 5*0(2), which are of the form 
t ^ for integer fc, \t\ = 1. When = 0, we have the trivial representation. Now since 
i?(l,2;6') is an element of 5*0(2) C SO{N), we can expression pa(i?(l, 2; 6^)) in terms of 
irreducible representations of 5*0(2), so for aj > 0, 

m 

Pa{R{l,2;e)) = J2(^jCOs{je). 

j=0 

Clearly, Ca(7r) = Y2j (yj{—iy > ao — «i. The a^-'s admit the following interpretation. By 
Weyl's character formula, one can express the restriction of pa from SO{N) to SO{N — 1) 
as follows: 

(2n+l) I n\ (2n) 

Pa '\sO{2n) = ^ Pb 

\bi\<ai<b2<...<b„<a„ 

J2n)| _ /T\ ^(2n-l) 

Pb \sO(2n-l) - Pa 

so each irreducible representation of 5*0 (A^ — 1) appears at most once when restricted from 
SO{N). This yields a combinatorial description of aj, namely the number of shrinking 
ensemble paths with the particle non-intersecting condition as given above such that the 
topmost particle arrives at ±j in — 2 steps (since we are going from SO{N) to 5*0(2)). 
We have the following two immediate consequences: 

Lemma 0.1. ai > a2 > ■ ■ ■, so the aj are monotone non-increasing starting at j = 1. 
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Proof. We will let /3j be the number of ensemble paths that leads to the af -* = j at the 
S0{3) level. Then condition on = j, j > 1, the next step one could get b\ E [—j,j], 
each with multiplicity 1. Thus if 1 < r < s G Z, the conditional multiplicity of a*^^-* = ±r is 
no less than that of a*^^-* = ±s. This shows > ag- Since there is only one copy of a^'^^ = 0, 
ao is not necessarily > ai. □ 

Lemma 0.2. Let he as defined in the previous proof. Then 2/3j — aj — Oij+i for j > 1, 
and Po = ao — ^ai. Furthermore we have /3j > 

Proof. The relation between ^ and a follows from the following two equations: 

m 

i=0 
m 



1=3 



These are consequences of the interlacing condition for the branching rule from SO{?>) to 
5*0 (2). Next we condition on a[^^ < ftg^'' < a^2^ . Then can be any integer in the range 
[—ai'\ ai^"*], each with multiplicity 1. So the number of branching paths down to S'0(3) with 

a\' — j is given by f3j := [1 + 2(j A a\ ')]I^^^j^(4)y where (3j :— ' ^ ' is the conditional 
/3 so to speak. Then it is clear by case analysis that 



0' > l±^/3' 



Now summing over all the conditions af^ < llf^ < a'^^ we obtain the second assertion. □ 

Observe now that X]j>o'^i ~ ^<^' since each shrinking ensemble path represents one copy 
of some irreducible representation of 5*0(2), which are all one- dimensional. We shall use the 
notation aj — aj/da, and similarly f3j = fij/da- Note that ^j>o/3j 7^ 1- 

The last lemma implies the following bound: 

Lemma 0.3. If ao > ai, then cto < O{{ao - ai)^/^) < C'(Ei>o "i(-lV)- 

Unfortunately, the condition in the above lemma only holds for small representations. 
Proof. Since aj — aj+i — 2/3j < 2(1 + 2j)/3o by telescoping, we have 

j 

ao — Qj+i = ao — ai + ai — a2 + . . . + aj — a^+i < 3^(1 + 2j)(3o, 
i.e., Oj+i > ao — {2j + 3)j(3o. Now using the fact that ^j>o <5j = 1, we have for all A; > 0, 

k 

{k + l)ao - Y,{2j + 3)j/3o = {k + l)ao - f{k)Po < 1, 



3 



where f{k) has cubic growth. So 



Optimizing over k on the right hand side, we obtain 

= c{ao - 2^1)^^^ 
< c{ao - aiY^^ 

Where in the last inequahty wc used monotonicity of aj, j > 1. □ 
Next we come back to bounding the total variation distance using the character sum: 



a^(0«) 



where (0") is the label for the trivial representation. We first give an upper bound on da- 
Since we know d(on) = 1, it suffices to look at the ratio: 

da ^ gg + g - 1/2 (g. + r - l/2f - (a, + s - l/2f 



dm t=i ^"^/^ r>s (r - 1/2)2 - (5 - 1/2)2 

Now suppose a = (0""*^, an-k+i, ■ ■ ■ , On), i-e., the first n — k coordinates are all 0. Then we 
have the following bound: 

J n I • 1 /o " I 111 

rfa _ -TT gj+J-l/2-rT -TT ttr + T — Sar + T + S — 1 

11 .-1/9 11 11 



d(Q■^^ j — ^l'^ T — s r + s — 1 

n n—k ^ n—k ^ 

< c TT exp(a^ ) exp(ar — ) 

r=n— fe+l s=l s=l 

<c TT ( ^a.^ r + n-k 

- 11 V-n + A;^ ^ r + 1 ^ 

r=n—k+l 

En 

< c 



1 Upper bound on the mixing time of the uniform 
Rosenthal walk 

Consider the following random walk on SO{n). Each step, an independent uniform angle 
^ G [0, 27r) is chosen and then an independent and uniform element in the conjugacy class of 
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R{1,2;6) is chosen. The walk then multiphes itself by this doubly stochastic element to reach 
the next step. This is closely related to the Kac random walk on SO{n), where instead of 
choosing a random element in the conjugacy class, one chooses a random pair of coordinates 
^ ^ i < j ^ n, and multiplies the chain by the random matrix R{i,j; 0). In the Kac case, no 
polynomial mixing time upper bound is known [8] (though see [9j, [1], and [3] for results in 
other modes of convergence). Here we present a polynomial mixing time argument for the 
corresponding Rosenthal walk. 
First we need a simple lemma 

Lemma 1.1. Consider a random walk K generated by fi on a compact group G which is 
constant on conjugacy classes and another one J generated by v which can be of any type 
(doesn't even need to be ergodic). Then the C? distance to stationarity satisfies the following 
censoring inequality: 

ll/i*^ * Z/ * /i*^ * . . . * Z/ * yU*'= — U\\2 < \\K^ — U\\2. 

where ti = t and ti > 0. 

Proof. This simply follows from the Plancherel's identity. It clearly suffices to show that for 
each irreducible representation p of G, 

Tr/i*i'*17*'/i*2(p) < Tr/I*(p). 

But this is true simply becuase fi{p) = cl for some constant c so can be collected in the front 
of the final matrix product. □ 

Now we will consider two auxilliary random walks to the uniform Rosenthal walk. Let K 
be the Rosenthal walk with 9 uniformly supported on [e, 2tt — e], and J be the unconditional 
walk, if we condition the Rosenthal walk on 6 in the range [e, 2tt — e], then we obtain K 
above. So the strategy now will be to show that K mixes in 0(n^log?2) steps and that with 
high probability in 0{n^ logn) steps there will be that many steps where K is taken. 

The second claim is easy to deal with. Let denote the number of steps in which 
e e [e,27r-e]. 



Pc*[#i^>t] = X:('/)(l-eye 



For c sufficiently large, the distribution of converges to a Gaussian centered at ct(l — 2e) 
with standard deviation about i/ct- So if we choose c > 2 then certainly > t with 
overwhelming probability. 

So now it suffices to bound the mixing time of K. 

• Regime 1. When r{n) = J^JLo^ji^^Y > V^, we can infer that > en? so the full 
monotonicity of the Fourier coefficients holds. Let 

T (,) .= (2^-^)' tH 

"^^^ ■ 22-1 (^^. + J _ 1/2) Ur^^dar + r - 1/2)2 _ (^^. + ^ _ 1/2)2) • 
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Prom the Rosenthal bound on r(7r), 

n 

we see that if a„ > 12n, then ra(7r) < 1/6. This is because the probabihty measure on 
[n] defined by 

:= Ton(j) - 



2^-1 (i - 1/2) n.^,[(r - 1/2)2 _ _ 1/2)2] ' 

is concentrated near j = 1 with a Gaussian decay of window size n^/^, as easily con- 
firmed by Stirling approximation. Furthermore for any a e (N")''", Ta{j) < ^x{j) using 
the monotonicity Oq < < . . . < a„. Therefore ^j>cni/2 2a0) = o(l), and we only 
need to focus on j < cn^l'^ . For the latter we simply compare Tj^f) with //(j) and 
bound the ratio uniformly: 



[a^ +j- 1/2)2 _ (^^^ ^ ^ _ 1/2)2 ^ ^ ^ j + ^ _ 1 



(j - 1/2)2 -(n- 1/2)2 - j+n-1 

So if a„ > 12n, this ratio is > 6. In fact this is true for all j G [n]. Therefore by 
monotonicity of expectation, we easily get rain) < 1/6. 

Now observe that a„ is precisely the number of nonzero a/s in the Fourier expansion 
of r{6), using the branching rule characterization. Therefore by monotonicity of a/s, 

Oil > — 

Q!l + Q!3 + • ■ ■ 



> 



> 



1 1 -r(7r) 
2 a„. 



Finally we have for cos^ < 1 — e (which is the case with our assumption that Q is 
bounded away from mod 27r), that 

r{Q) < 1 - (1 - cos^)ai < 1 - (1 - cos^)^^— 

2a„ 

Now if we let A be all the irreducible representation indices of SO{2n + 1) with the 
property that ra(7r) > 1/6, then Rosenthal clearly implies that 

5^4ra(7r)"i°^" = o(l), 

which in turn gives 

5^4raW™"'°^" = o(l). 

agA 
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Regime 2. From Rosenthal's work [TO], we know the character ratio at a particular 9 
is bounded as follows: 

{2n-l)\ ^ 1 



(2sin^/2)2-i ^ a,- 1 H.^.K + r)2 - {a, + j] 



Let us denote Wg = Y[r>s I ''"'^^''rZ-l""^'^'' I • Using the fact that the character ratio is 1 
for the trivial representation, taking the smallest Wq gives a bound of r{6) as follows: 

n 

r(e) < y — . 

^ ' - (sin^/2)2«-i ^ 

We also have the control over the gap between the smallest and largest Wg, max^ Ws < 
(miug Ws)'^"'. This follows from the simple and crude observation that for any fixed j, 

1 

2' 

(a„ + aj + n + j) > (a,. - + r - s). 

Actually the bound above is the best one could hope for, as it is achieved by the 
following family of representations (0"~^, k), as k ^ oo. 

But now we have r{9) < (sipg/2)2"-i Wj)~^n, whereas the dimension d < (maxj ly^)^"/^ 
by throwing away the corresponding terms for the trivial representation, and estimat- 
ing the factor Hgl'^g + 9 ~ 1/2) trivially by the product of the W/s. So denoting by 
Wo the minimizer of Wj, we need t such that the following expression is small: 



{ttn + ttj + n + j) > -{ttr + tts + r + s) 



It is important that our upper bound for r{6) is < 1. to accomplish that we need 
Wq > e^^'^\ In that event, we certainly can take t = ^l{n'^) and get the required 
bound. 

Regime 3. Now we need to consider the case when our upper bound for r{6) > 1. 
One possibility is taken care of by Regime 1. So it suffices to consider the situation 
where r(7r) < 1/6. But that means < ^ by an easy computation using the fact 
that ao — tti + ^2 < r(7r). Integrating over 9 £ [e, 2tt — e], we obtain 



2n - 2e 



/ ri9)d9 = ao-J2- rsin(ej) 



It is clear that < 

So for fixed e G (0, vr — 1), we have the uniform bound Jj^ 2TT-t] ''^{(^)d9 < 1 — 0^(1 — «o)- 
So in other words, r{9) < 1 — 5 for some 6 > 0. On the other hand, the dimension is 
bounded by 

So taking t = Q{n^) suffices to kill all the harmonics in this regime. 
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In summary, the mixing time of the uniform Rosenthal walk is 0{n^). There is a cor- 
responding lower bound in £^ of order (see Porod and Rosenthal). The above argument 
clearly shows that the mixing time of the uniform walk is of the same order as for a fixed 6 
(or uniform ^ is a region bounded away from 0). So the conjectured total variation mixing 
time is still nlogn. 
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